ON THE CONVERGENCE OF AXIALLY SYMMETRIC VOLUME PRESERVING 

MEAN CURVATURE FLOW 



MARIA ATHANASSENAS AND SEVVANDI KANDANAARACHCHI 



ABSTRACT. We study the convergence of an axially symmetric hypersurface evolving by volume- 
preserving mean curvature flow. Assuming the surface is not pinching off along the axis at any time 
during the flow, and without any additional conditions, as for example on the curvature, we prove 
that it converges to a hemisphere, when the hypersurface has a free boundary and satisfies Neumann 
boundary data, and to a sphere when it is compact without boundary. 



1. Introduction 

Consider n-dimensional hypersurfaces M t , defined by a one parameter family of smooth immer- 
sions x t : M n — > M. n+1 . The hypersurfaces M t are said to move by mean curvature, if x t = x(-,i) 
satisfies 

(1.1) ^-x(p,t)=-H(p,t)v(p,t), p£M",t>0. 

at 

By v(p,t) we denote a designated outer unit normal of M t at x(p, t) (outer normal in case of 
compact surfaces without boundary), and by H(p, t) the mean curvature with respect to this normal. 

If the evolving compact surfaces M t are assumed to enclose a prescribed volume V the evolution 
equation changes as follows: 

(1.2) ^- X (p,t) = -(H(p,t)-h(t))v(p,t), P GM n ,t>0, 
at 

where h(t) is the average of the mean curvature, 

hh d 9t 

and g t denotes the metric on M t . The surface area \M t \ of the hypersurface is known to decrease 
under the flow (see [2]). 



We are interested in an axially symmetric surface, which encloses the volume V, and which has 
a nonempty boundary contained in a plane II that is perpendicular to the axis of rotation. Motivated 
by the fact that the stationary solution to the associated Euler Lagrange equation satisfies a Neumann 
boundary condition, we also assume the surface to meet that plane II at right angles along its 
boundary. Assuming the surface to be smooth, it will also intersect orthogonally the axis of rotation. 

We consider the case where the surface will not pinch-off along the axis of rotation during the 
flow, having only one intersection with that axis at the point that is the furthest from the supporting 
plane IT, and prove that the surface converges to a halfsphere. 
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The methods we use also apply in the case of an axially symmetric surface without boundary 
having a similar lower height bound, and in that case we prove in Section 8 that the flow converges 
to a sphere. 

The results in this paper make use of the axial symmetry, and no additional conditions on the 
curvature of the surface are assumed. Converge to spheres has been previously proved for the 
volume flow by Huisken in [7], for compact, uniformly convex initial surfaces; and in [10], Li 
assumes bounds on the traceless second fundamental form. 

Our results can be seen as complementing the work of the first author in [2, 3], and the second 
author's work on her PhD dissertation: In the case of the surface behaving like a "bridge" between 
two parallel surfaces, if one were able to flow through singularities, the axially symmetric volume- 
preserving flow would converge to a number of spheres and (possibly) two hemispheres on the 
parallel planes, like beads strung along the axis of rotation. 



2. Notation, definitions and assumptions 

In the case of the surface M t intersecting the obstacle II, we will at different stages divide it into 
two parts as in [1]: one adjacent to the plane and the remainder that contains the (only) intersection 
with the axis of rotation. 

Let II = {(xi, . . . ,x n+ \) G IR n+1 : x\ = 0} and M t be contained in the right halfspace, 
M t C {x\ > 0}. We use R t as the generic notation for the part of the surface closest to the plane, 
and Ct for the rest - the cap that intersects the axis of rotation - , and we will introduce various 
superscripts depending on the situation that will be made clear in the text. 

We denote by P(t) = (d(t),0) the "pole": the point of intersection of M t with the axis of 
rotation. We assume that there are no singularities developing, so that P(t) is the only point of 
intersection of M t with the axis of rotation for all time. We are interested in those solutions where 
the generating curve of the initial hypersurface is smooth and can be written as a graph over the xi 
axis except at the pole. 

We use the notation 

Pt : [0, d(t)] -> R 
for the radius function of the surface of revolution. 

Let ii, ■ ■ ■ , i n +i be the standard basis in W l+1 and ii be the direction of the axis of rotation. We 
denote the quantities associated with the cap with a tilde ~ , and in this context we work with the 
vertical graph equation. 

Furthermore we define the following quantities on M t : 
Let lj = € IR n+1 , x = (0, X2, ■ ■ ■ , x n+ \), denote the outer unit normal to the cylinder intersect- 
ing M t at the point x(p, t). We call u = (x, oS) the height junction of M t , and set v = (v,^)^ 1 . 
Note that v corresponds to y/l + p 2 , and will be used to obtain gradient estimates. 
The respective quantities on the cap Ct are the height measured from the plane II, u = (x, ii), and 

v = (^hr 1 . 



ON THE CONVERGENCE OF AXIALLY SYMMETRIC VOLUME PRESERVING MEAN CURVATURE FLOW 



3 



We cut the hypersurface in two regions using the plane L a (t), which is parallel to IT where (u, ii ) \ La (t)nM t = 
1/a , with a being a constant. We define the cap, determined by the inclination angle, as the con- 
nected component of M t containing the pole P 

C t a =|x(p,*)€M t :i<(i/,ii)<l| , 

and we call Rf = M t \Cf the cylindrical part of the surface. Note that L a (t) is chosen such that 
the specific inclination angle is achieved nowhere else between that plane and the pole P(t). As 
long as the flow is smooth, Cf is by definition a graph over the X\ axis except at the pole. 

Assumption 2.1. We assume for any a > 1 there exists a constant c(a) > depending only on a 
such that u \r<* > c(a) , i.e we assume a lower height bound in Rf , independent of time, dependent 
on a . 

Thus P(f) is the only point of intersection of M t with the axis of rotation for all time. The assump- 
tion prevents singularities developing on the axis of rotation. 

For an axially symmetric surface the mean curvature is given by 

h= — p -^ + ^K, 

(l + /9 2 )2 p(l + p 2 )2 

while the principal curvatures are k = - — *- and p = — } We also introduce another 

quantity q = (is, ii) u' 1 such that in particular p 2 + q 2 = u~ 2 . 

3. Height Estimates 

In this section we prove that M t satisfies uniform height bounds: both, the height function u 
as defined above, and also the height when measured as distance from the obstacle IT, i.e. u, are 
bounded. 

Lemma 3.1. The evolving surfaces M t satisfy the uniform height bound 



u < R = 




Proof. We follow a method as in [2] in getting bounds for u. Assume there exists an R such 
that UM t > R at some given time t. Since the surface area is decreasing under the flow, and by 
comparing to the projection of the surface onto the plane, we have 

\M \ > \M t \>uj n R n , 

where uj n is the volume of the n dimensional unit ball. Therefore 




would contradict the fact that the evolution decreases the surface area. □ 

Lemma 3.2. There is a constant I such that the evolving surfaces M t satisfy the height bound 

u < I , 
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that is the distance from the plane IT is uniformly bounded. 




Figure 1. The cylinder of radius c(a). 

Proof Here a = and 

C t a = |x(p,t)€M t :^<( I /,ii)<l| , 

and Rf = M t \Cf. From the Assumption 2.1, we know that u > c(a) in Rf . As u \dc a < R an d 
\f/\ > tan (f - 9) = in C? we have 

d{t)-u \ dC f< Rtan6 = R^/a 2 - 1 . 

Assume there exists a length l\ such that u \r<*> l\. Then 

|M | > \M t \ >nuj n c n - 1 {a)h, 

where now we compared \M t \ to the surface area of an n dimensional cylinder of radius c(a) and 
length l\. Having l\ > s would contradict the fact that the evolution decreases the surface 

b 1 fc 1 nw n c n 1 (a) 

area. Therefore 

u< . + Ry/^l=:l. 

nu n c n 1 (a) 

□ 

Next we show that the length of the generating curve is bounded. 

Lemma 3.3. Assume Mt to be a smooth, rotationally symmetric hypersurface, with a radius func- 
tion p(xi,t) > Ofor x\ € [0, d(t)). Then there exists a constant c*, such that 

rd(t) 



/•' . • ■ 

/ v 1 + p' 2 dx\ < c* , 
Jo 



independent of time. 



Proof. Let us divide M t into Rf and Cf for any a > 1. As the surface area is decreasing under the 
flow 

|M t | < |M |, 

rd(t) 



/•aw 

2vr / ^Vl + p^xi < |M |, 

■/ 

2vr / p^Vl + P' 2 dxi < 2tt / P n ~ 1 y/l + f/ 2 dx! < \M \ . 
Jo Jo 
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From the Assumption 2.1 

rL a (t) 



2irc n - 1 (a) / V 1 + p' 2 dxi < | M 1 , 
Jo 

fL a (t) 

/ y/1 + p' 2 d Xl < 

JO 



\M \ 



Iq 27rc n - 1 (a) ' 

We can estimate the length of the generating curve of the cap Cf by I + R. Therefore 

rd(t) | Mq | 



/ y/l + p' 2 dxi < 



27r C «-i(a) +/ + jR " :C *' 

□ 



We now derive an a priori estimate for h(t) for any solution of the graphical equation. 



4. Estimates on h 

Lemma 4.1. Assume M t to be a smooth, rotationally symmetric hypersurface, with a radius func- 
tion p(xi,t) > Ofor x\ G [0, d{t)). Then there is a constant c\ such that < h(t) < c\ throughout 
the flow. 

Proof. Following [3] we parametrize Mt by its radius function p G C7°°([0, d(t))), then clearly 

H = P -^ + 



(1 + p 2 )^ p(l + p 2 )5 

From Lemma 3.3, we know that \Jl + p 2 dx\ < c*. Our proof follows the ideas of [2] the 
difference being the boundary term when integrating by parts. For the sake of completeness we will 
include it here. For the second term of 

h ® = Tin f ( k + ( n ~ l )P) d 9u t € [0, T) , 

\ M t\ JM t 

we have 

n-1 f d ^ „_ 2 {n-l)R n ' 2 l 



o<^^y o p n ~ 2 (x 1 ,t)d Xl < 



n\ jo \M t 
since p < R and d(t) < I by Lemma 3.1 and 3.2. 

p - d 



For the first term note that jj+^rj = (arctan p). Therefore 

f f d ^ d 

\ kdgt = — — — (arctan p)p n dx\ 
JM t Jo dx l 

rd(t) 

(4.1) = (arctan p)p n 1 \ Xl =o —(arctan p)p n 1 \ xi=d n\ +{n — 1) / (arctan p)pp n 2 dx\ 

Jo 

rd(t) 

= (n — 1) / (arctan p)pp n dx\ , 
Jo 
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as arctanp = when x\ = 0, and p(d(t)) = at the pole. As < (arctanp)p < ||p| < 

-1)J 

\M t \ 



I y'l + p 2 we obtain 



2 Jo 

(n- l)^^"- 2 ^ 

- \m\ 2 ' 

where we have used Lemma 3.3. 

From the isoperimetric inequality and the fact that the flow decreases surface area we know that 

n 

V"+i < c \M t \ < c\M Q \. 

Hence combining these arguments we conclude 

[Hdg 

< ^r-^ < Cl. 

□ 

5. Evolution equations and gradient estimates 

The maximum principle for non-cylindrical or time dependent domains is discussed in [9]. We 
use that version of the maximum principle in this paper. 

Lemma 5.1. We have the following evolution equations: 



« (1 
© (1 


- A)u = 


h n—1 . 

v u ' 


-A)u = 


k ■ 


(iii) (1 


-A)v = 


-\A\ 2 v + (n-l)%-l\Vv\ 2 ; 


(iv) (| 


-A)v = 


-\A\ 2 v - i\Vv\ 2 ; 


W a 


- A)H = 


-(H-h)\A\ 2 ; 


(vi) (I 


-A)|^| 2 


= -2\VA\ 2 + 2\A\ 4 -2hC; 


(viO (| 


-A)p = 


\A\ 2 p + 2q 2 (k - p) - hp 2 ; 


W (| 


- A) k = 


\A\ 2 k - 2(n - l)q 2 (k - p) - hk 2 



where C = g t] g kl g mn hikhi m h n j, with g l i denoting the components of the inverse of the first fun- 
damental form, and hij those of the second fundamental form. 

Proof, (i) and (iii) are proved in [2], (v) and (vi) are in [7]. 
(ii) For u = (x, ii) we have 

Tt* = (jt*> h ) = - {H - h){u > hh 

and 

Au = (Ax,ii) = 

so that 

± - A) u = h(v,h). 
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(iv) For v = (is, we have 

The evolution equation follows from the following well known identity (see for example [5]) 

Av = -i 2 (VH, ii) + v\A\ 2 + 2v~ 1 Vv 2 . 
Following the same way as in [8], for (vii) we have , 

d d . _ 2 2\l/2 

dt p= dt iu - q) 

= Ap + p^lVpl 2 + p^\Vq\ 2 - 3p _1 u _4 |Vu| 2 + p~ 1 u~ 4 
- qP' 1 (\A\ 2 q + q(p 2 ~q 2 - 2kp)) - hu~ 2 + hq 2 . 
Equation (vii) follows from the following relations: 

(5.1) ViU = 5nqu, Vi {v, ii) = kpu , V \q = 5a(q 2 + kp) , 

(5.2) V i p = S il q(p-k), \A\ 2 = k 2 + (n - l)p 2 , u~ A = p 4 + 2p 2 q 2 + q 4 . 

The evolution equation for H was derived in [7] , and (viii) follows from (v), (vii), and the fact that 
H = k + (n — l)p. □ 

We proceed to obtain gradient estimates in the different parts of the surface: for the cap using 
the vertical graph equation and part (iv) from above Lemma 5.1, and for the cylindrical part away 
from the cap, using the evolution equation (hi) in Lemma 5.1. 
The quantities u and v are used on the cap. 

Lemma 5.2. The gradient estimate 

v < a 

holds on the cap Cf. In addition there is a constant C2 (a), such that 

v < c 2 (a) 

for the cylindrical part Rf. 
Proof. Note that 

so that by the maximum principle v < max(maxq v, m&XQc? v). By definition in Cf we have 
v < a, and this is supported by the evolution equation! 

From the assumption we know that u > c(a) in Rf. As in ([2], Proposition 4) we calculate, 

- A^j u 2 v = -\A\ 2 u 2 v + (n - l)v + 2uh - 2(n - l)v - 2v\Vu\ 2 

2 _ 

VvV(u 2 v) 

v 

< 2hu - (n - l)v . 
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If v > the right hand side is negative, and proceeding as in [2] we conclude v < 02(a) in Rf. 
It is important to note that on the boundary of Rf, either v = 1 (along the intersection with II), or 
v = ^/=§=p Thereby, we have bounds for v and v in Rf and Cf respectively. □ 

Remark 5.3. These gradient bounds guarantee that Rf remains a graph. As Rf remains a graph 
for all a > 1, we see that M t \P(t) remains a graph as well. 

As the height of the graph is bounded, for the minimum d(t) we find a lower bound from 

fd(t) j-d(t) 
V = / uj n p n (x)dx l < uj n R n / dxi = uj n R n d(t) . 
Jo Jo 

Lemma 5.4. Let xo(i) be a boundary point ofC^, which without loss of generality we can assume 
to lie on the generating curve and such that (i/(xo(i)), ii) = ^ (with some abuse of notation for 
the corresponding normal v(xo(t))). Then H(xo(t)) > OforO <t< T rnax < 00. 




Figure 2. x (i) - the boundary point of cf 1 . 

Proof. Suppose H(xo(t)) < 0, then by continuity there is a connected region C^ H C cf 1 , with 

xo(t) G dCt , which clearly can be chosen to be axially symmetric, and such that H | ^ h- < ~- 

c t 

0. Let xi (t) denote the other boundary point along the generating curve in c/^' H C c/^, and let 
a(t) = (xo(t), ii), b(t) = (xi(t), ii) denote the x\ coordinate of xo(t), xi(t), respectively. Then 



Icf' H ~ 



b(t) 



P_ „n-l 1 /„ i\ „n-2 

la(t) V l + P 2 

The second term being positive, that means that the first is negative, and given the bounds on the 
radius we find 

/ I k dx\ = I I — : — (arctan p) dx\ < . 

Ja(t) V 1 + PV Ja(t) V / 

This results in arctan p(a(t)) < arctan p(b(t)) and by the choice of a(t), —j < arctan p(b(t)). 

/2 

However, this is not possible in C t , where — f < arctan p < — ^, contradicting our assumption 
and therefore H(x (t)) > 0. □ 



> / K _ Hdg = I ( -r^P^ + (n - W ) d Xl 
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6. Curvature estimates 

Proposition 6.1. There is a constant C2 depending only on the initial hypersurface, such that | < 
C2, independent of time. 

Proof. We calculate from Lemma 5.1 

d ( k\ A k 2„ „ fk\ a 2 . ., . ,, hk , 

- - =A- + -V iP V i - +2% (p-fc n-lp + fc + — (p-fc • 

at \p ) p p \P J P V 

If | > 1 then M (p _ jfe) < o. This implies that 

k ( k 

(6.1) — < max l,max — 

Note that for this consideration, the smooth function k/pis defined over the whole surface, and in 
view of the orthogonality on the boundary, via a reflection argument there are no boundary data 
involved. □ 

Proposition 6.2. There exists a constant C3 such that 

\A\ 2 < cs. 

Proof. We proceed as in [6] and [2] and calculate the evolution equation for the product g = 
\A\ 2 ip{v 2 ) in Rf 2 , where ip(r) = A _^ r , with some constants A, fi > and v = (v,oj}~ 1 . From 
the evolution equation of g we find the inequality 

j t ~ A) g < -2,,g 2 -2X^Vv.Vg- {x l X ^ 2)2 \Vv\ 2 . 

We estimate the second last term as in [2] using Young's inequality and obtain 

-2hOp(v 2 ) < 2h\A\ 3 <p(v 2 ) 

<§WW) + ^VV) 
We choose fi > | and A > /imaxw 2 . As ip'v 2 = pp^rpV? we have 

As u > c = c o in we get 

2(n-l)A 



u 2 {\ — /IV 2 ) 

Therefore we have 

d . A 2 



9 <ca9- 



; ^ - Aj g < -c 5 g + c 6 g - c 7 Vi; • V# + c 8 (/i, max v) 
<-c 5 (g--^-j - c 7 Vv -Vg + c 9 . 
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When g > ^ + we get the right hand side to be negative. On dR/ 2 we have H = k + (n — 
l)p > by Lemma 5.4 and also as | < C2 we get < c on the boundary point. Here we have 

\A\ 2 = k 2 + (n - l)p 2 < (c 2 + n - l)p 2 < C/>~ 2 < Ccq 2 . 
From the maximum principle 

g < max max g, max |74| 2 (/?(f 2 ) 

As 2 ) is bounded as v is bounded we have a bound for g in r/ 2 - If we calculate the evolution 
equation for g = \A\ 2 cp(v 2 ) on c/ 2 then we get the same evolution equation without the last term 
on the righthand side. Thereby we get a bound for g in the same way as above. □ 

Proposition 6.3. For each m > 1 there is C m such that 

\V m A\ 2 <C m , 

uniformly on M t for < t < T max < oo. 

Proof. Having obtained uniform bounds on \A\ 2 and h the proof is a repetition of that of Theorem 
4.1 in [7]. □ 

Thus we have long-time existence for the flow. 

Corollary 6.4. 

7. Convergence to surfaces of constant mean curvature 

This is as in [2] : Having long-time existence, Proposition 8 of [2] gives convergence to a con- 
stant mean curvature surface, which in our case is axially symmetric. By the classification of the 
Delaunay surfaces [4] it has to be a half-sphere. 

8. Other convergence results 

Using the same estimates with very few changes one can show that a compact, axially symmetric 
surface without boundary, which encloses a volume V and intersects the axis only at two endpoints 
throughout the flow, will converge to a sphere. We will only explain the parts that are different from 
the previous result. 

8.1. Height Estimates. The height estimates 3.1 and 3.2 change as follows. 
Lemma 8.1. The height function u satisfies u < R = (^j^J ". 

Proof. Assume there exists an R such that uu t ^ R at some given time t. Take a plane perpendic- 
ular to the xi-axis and intersecting the surface. This plane divides the surface into two parts, and 
by projecting both parts onto the plane we find 



|M | > \M t \ > 2oj n R n . 
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Taking 



2u n 

would contradict the fact that the evolution decreases the surface area. □ 

Lemma 8.2. Let e(t) and d{t) denote the two tips of the surface on the left hand and right hand 
side respectively. Then d(t) — e(t) < I = ni ^ M ^-i + 2it 

Proof. As in Lemma 3.2 let a = t^tts. From the Assumption 2.1 we know that u > c(a) in Rf . 
As u \ dc a,i< R and \p\ > tan (| — d) in C"' 1 for i = {1, 2} , we have 

d(t) — u \g C a,i< Rt&n9 = R\f a 2 — 1 , and 



u Ig^a.a — e(t) < Rt&nO = Ry a 2 — 1 . 

Assume there exists a length l\ such that u \ Rf > l\. Then by the previous argument again we can 
say 

|M | > \M t \ > nuj n c a ~ l (a)h , 



where now we compared \M t \ to the surface area of an n dimensional cylinder of radius c(a) and 
length h. If h ; 
area. Therefore 



length l\. If l\ > s this would contradict the fact that the evolution decreases the surface 

& 1 1 nu n c n 1 (a) 



n<^M T + 2i? V / ^T. 
nuj n c Q 

□ 

Lemma 8.3. Assume Mt to be a smooth, rotationally symmetric hypersurface, with a radius func- 
tion p(x±,t) > Ofor x\ € (e(t), d(t)). Then there exists a constant c*, such that 

rd(t) 



ra x , 

/ a/1 + p' 2 dx\ < c* , 
Jo 



independent of time. 



Proof. The proof is the same as in Lemma 3.3 after taking into account the two caps on either side. 
Here we would have 



d(t) | Mq | 



□ 



Lemma 8.4. (Estimates on h) Assume M t to be a smooth, rotationally symmetric hypersurface, 
with a radius function p{x\,t) > Ofor X\ G (e(t),d(t)). Then there is a constant c\ such that 
< h(t) < C\ throughout the flow. 

Proof. The only change to the proof of Lemma 4.1 would be in the boundary values 4.1. Here the 
new boundary values would be 

(arctanp)//^ 1 \ Xl=a{t ) -(arctanp)/?™- 1 \ Xl=b{t ) 
As p(a(t)) = p(b(t)) = 0, the boundary terms dissapear and we get the same estimate for h. □ 
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Lemma 8.5. (Gradient estimates) The gradient estimate 

\v\ < a 

holds on the caps C"' 1 , i = 1,2. In addition there is a constant c, such that 

v < c 

for the cylindrical part Rf. 

Proof. The gradient estimates are as in Lemma 5.2, but in this setting instead of one cap Cf we 
have two caps on either side, and the same estimate would suffice for both caps. □ 

Concluding this section, we remark that H > at points where the caps cf^' 1 , i = 1,2, 
meet the cylindrical part Rf* of the surface. The proof is using the same arguments as the one 
for Lemma 5.4 after the appropriate adjustments of the sign of arctan p for cap on the left of the 
surface. The results on curvature estimates and the convergence to a limiting surface of constant 
mean curvature follow along the same lines as previously proved. In this case the limit surface is a 
sphere. 
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